NSF-ITP-98-089 
CERN-TH/98-298 
hep-th/9809106 



Aspects of QCD Dynamics from String Theory 

Akikazu Hashimoto and Yaron Oz 6 

"Institute for Theoretical Physics 
University of California 
Santa Barbara, CA 93106 

fe Theory Division, CERN 
CH-1211, Geneva 23, Switzerland 



Abstract 

We study dynamical aspects of large N QCD4 via supergravity on Anti de Sitter 
black hole geometry. We compute the mass spectrum, the topological susceptibility 
and the gluon condensate and make a comparison to lattice simulations. We present 
analogous results for QCD3 . 
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1 Introduction 



Maldacena [1] conjectured the existence of a duality between large N superconformal 
field theories and superstring or M theory on Anti de Sitter (AdS) backgrounds with R-R 
fluxes. This conjecture reveals interesting properties of the superconformal field theories, 
in particular in the strong coupling region. In [2] Witten proposed an extension of the 
duality to non-supersymmetric theories such as pure QCD. The idea was to heat the 
superconformal field theories and break the superconformal invariance. The AdS space is 
then replaced by the Schwarzschild geometry describing a black hole in the AdS space. 

When the curvature of the space is small compared to the string scale and the Planck 
scale, supergravity provides an adequate effective description of string and M theory 
respectively. This description has been extensively studied recently. It exhibits qualitative 
agreements with pure QCD in three and four dimensions [3]. In particular confinement 
and mass gap are two of the important properties of the supergravity description [2]. This 
leads to the hope that we are at least describing a theory which is in the same universality 
class as QCD. 

One can also perform quantitative computations via the supergravity description. For 
instance, the spectrum of the scalar glueball masses has been computed resulting in a 
surprising agreement with lattice simulations [4]. The supergravity limit of string theory 
is the limit where the string tension is infinite. This corresponds to the strong coupling 
limit of the gauge theory. The strong coupling expansion in the gauge theory is the a' 
expansion of string theory. It is therefore clear that a full solution of QCD requires a 
detailed knowledge of this expansion, which means an understanding of string theory 
with R-R fields on these backgrounds. Since supergravity contributes the leading term in 
this expansion it is important to have a detailed understanding of it. We will take this 
viewpoint in this paper and use the supergravity description to study various aspects of 
QCD. Note, however, that an important assumption that is made is that there is no phase 
transition as a function of the 't Hooft parameter g\ M N . 

The paper is organized as follows: 

In section 2, we describe in detail the supergravity background and the conjectured cor- 
respondence to QCD 4 . 

In section 3 we compute glueball mass spectrum. 

In section 4 we compute the topological susceptibility. 

In section 5 we compute the gluon condensate. 
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In section 6 we make a detailed comparison with lattice and field theory results. 

Section 7 is devoted to a discussion and conclusions where we also present analogous 
results for QCD 3 . 



2 QCD Via String Theory 

The starting point for the construction of QCD 4 is the superconformal theory in six 
dimensions realized on N parallel coinciding M5-branes. Compactification of the theory on 
a circle of radius R\ gives a five- dimensional theory whose low-energy effective theory is the 
maximally supersymmetric SU (N) gauge theory with gauge coupling constant g\ = 2nRi. 
In order to describe QCD 4 , one compactifies this theory further on another S 1 of radius R 2 . 
The gauge coupling constant gyu in four dimensions is given by g\ M = g\jRi = 2nRi/ R 2 . 
In order to break supersymmetry, one imposes the anti-periodic boundary condition on 
the fermions around the second S 1 . 

As conjectured in [1], the large N limit of the six-dimensional theory is M theory on 
AdS 7 x S 4 . Upon compactification on S 1 x S 1 and imposing the anti-periodic boundary 
conditions around the second S 1 , we find M theory to be on the black hole geometry 
[2]. Taking the large N limit while keeping the 't Hooft parameter 2ir\ = gy M N finite 
requires taking Ri <C R 2 . In this limit, M theory reduces to type II A string theory and 
the M5 brane wrapping on S 1 of radius R\ becomes a D4 brane, and the large N limit of 
QCD 4 becomes string theory on a black hole geometry. 



2.1 The Supergravity Background 

The action for the eleven dimensional supergravity is 1 



S = 



The reduction of (2.1) to ten dimensions on a circle of period gives 



where 



K \ Q = (N/X)Tk 



ii- 



(2.1) 



(2.2) 



(2.3) 



1 Thc N = 1 supergravity action in eleven dimensions also contains Chern-Simons terms and a bound- 
ary extrinsic curvature terms which we omitted. 
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The near horizon limit of the black M5-brane metric is given by [5, 6] 



7 2 y 
ds = - 



H = 3V6L 3 dQ 4 



where 



L — N 



i=i 



3 K ii 



yo 



4ttL 2 T 



(2.4) 
(2.5) 

(2.6) 



2 7 7T 5 ' ' 3 

and dVL$ is the volume form on the 4-sphere. Compactifying x 5 on a circle and using the 
M theory/Type IIA duality one gets the Type IIA background corresponding to a non 
extremal D4 brane with worldvolume coordinates x 1 , ..,x^,t 



y° 



(i-|)^ 2 + E^ 2 
y i=i 

where 6 is the dilaton field. 



±L,_ yt 



+ —(l-^)-Uy 2 + Lydnl 



if, P.7) 



2.2 Fluctuations 

As will be discussed in the next subsection, in order to compute the correlation functions 
of QCD operators using super gravity /gauge theory correspondence, we need to examine 
the dynamics of small fluctuations of the supergravity fields. We will be interested in 
scalar operators of QCD4. The relevant fields of Type IIA supergravity are the dilaton 0, 
the volume factor, and the time component of the RR 1-form Cq. There will be a mixing 
between the type IIA dilaton and the volume factor which must be disentangled. A useful 
parameterization of fluctuations f(xi,y), g(x±,y), and h(xi,y), which in the language of 
eleven dimensional supergravity is appropriately diagonal, is 



ds 2 = 



4 
5"' 



\ 1 / 2 2 \ 4 

9(xi, y) Y, gudx l dx l + 1 - 7/(^1, y) - -g(xi, y) Y, gadx % dx % 

+ (l + ^f(x u y) - ^g{x u y)^j g 55 (dx 5 + Cdx^ 



1 - -g(x u y) g yy dy + (1 + g(x u y))g nn dVt 



(2.8) 



where g is the background metric. We consider the case C = h(xi,y) and C^o = 0. 
Without loss of generality, f(xi,y), g(xi,y), and h(xi,y) are assumed to depend only on 
Xi and y. In the language of Type IIA supergravity, these fluctuations are parameterized 
according to 



ds 2 



3 6 \^ { 1 6 \ 

1 + ^f(xi,y) - ^g(xi, y)j Y 9udx l dx t + [1 + -f(x u y) - -g(x u y)j Y gudx'dx 1 
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3 6 \ f 3 3 \ 

i + g/fai, y) - -^g(xu y)J g yy dy 2 + ( 1 + ^ /(^ v) + ^g( x u y)J gnndtt 2 

e 4, = (l + ^/(zi,v)-|/(zi,v))e* 

C = h(x u y). (2.9) 
By directly substituting these fluctuations into the supergravity equations of motion, 

iV - V* + V = 0, (2.10) 



we see that they satisfy the equations (r — y/yo) 

9M 2 

/ Qfl/f2 



r 4 - r~ 2 )f" (r) + (7r 3 - r" 3 )/'(r) + ^^/(r) = 0, 

/ 9M 2 \ 
(r 4 - r" V(0 + (7r 3 - r-V(r) + [^^g(r) ~ 72r ") = °> 

(r 4 - r- 2 )/i"(r) + (7r 3 - 7r->'(r) + f^h(r) = 0. (2.11) 
2.3 Field/Operator Correspondence 

The string/gauge theory correspondence asserts that the string modes (p couple at the 
boundary to gauge invariant operators O of the gauge theory. The generating functional 
for the correlation function of O is the string partition function evaluated with specified 
boundary values (po of the string fields. When the supergravity description is applicable 
we have 

( e -/**M*)0(*)) = e - IsG (ip ), (2.12) 

where Isg is the supergravity action. 

The supergravity fluctuations /, g, h couple to operators on the D-brane world volume 
in flat background. For small 't Hooft parameter, this coupling can be inferred from the 
form of the Born-Infeld action. These same operators, in the strong coupling limit, are 
described as "sources at infinity" in the dual AdS description. 

The Born-Infeld Action describing the D4 brane worldvolume theory is given by : 
Ibi = TrT 4 J d 5 x (e~^ y/ det(G ^ + + G^d.X^Xj) + V^C^^) , (2.13) 
where % — 5 . . . 9 and T 4 is the D4 brane tension. 

Dimensionally reducing along the t direction on a circle with period 1/T leads to 

Ibi = Tr^ j d 4 x (e~* '^^(G^ + + + \e^C Q T^T c ^ , 

(2.14) 
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where now, i = 0, 5 ... 9. In order to restore the canonical dimensions of Yang-Mills 
fields, we rescale F and X: 



^ = T~ 1 F, X = T~ 1 $, (2.15) 

where T' 1 = 2nl 2 . 

Expanding the action (2.14) to quadratic order gives 

I = T±f d'x (^e^ V ^ V / d^> 2 + _^_^C F^) , (2.16) 

where we dropped the scalars since they are expected to decouple from the dynamics in 
the far infra-red. Using [7] 

T p = ( ? - 1 (27r)( 1 -^/ 2 T( p+1 )/ 2 , (2.17) 

we get 

% 1 " (2.18) 



T 2 T (2n) 2 gl s T 2vrA' 
where in the last equality, we set 2ng s l s = (A/iV^T -1 which follows from M theory/Type 
IIA duality. Before substituting these fluctuations into (2.16), note that the metric fluctu- 
ation due to field f(x 1: y) is polarized along (3(dx 1 ) 2 + X^ =2 (cfc*) 2 ) which is not isotropic 
in the 1234 plane. One could however think of this as contribution from two pieces 
(J2i = i{dx 1 ) 2 ) and (dx 1 ) 2 . The second term must then couple to Tii(pi) component of the 
stress-energy tensor which vanishes due to conservation law p^T^ = 0. We will therefore 
think of fluctuation f(xi,y) as giving rise to a fluctuation of the metric of the form 

ds 2 = (l + ^f(x!,y) - -g(x!,y)^ YLgudx l dx l 

along the 1234 plane. Substituting these fluctuations into (2.16), we find that f(xi,y), 
g(xi,y), and h(xi,y) couple to the world volume operators via 

1 = Tr / A (i (£x) ( l ~ l f ^ y) ) p2 + i (£x) k™^***) ■ < 2 - 19 ) 

Note that all dependence on g(x 1: y) canceled out to this order. Indeed, the AdS/SCFT 
correspondence in six dimensions imply in view of (2.11) that g(xi,y) couples in six 
dimensions to an operator of dimension 12 in the (0,2) theory. This operator is where 
H = dB, which reduces to a dimension 8 operator in the four dimensional theory in the 
low energy limit. Comparing (2.19) to the Yang-Mills action 



'™^J d ' x {ik F2 + ^ Fp )' (2 ' 20) 
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jPC 


1 


2 


3 


4 


5 


6 


f(xi,y) 


0++ 


9.85 


15.6 


21.2 


26.7 


32.2 


37.7 


g(xi,y) 


0++ 


22.5 


28.8 


34.9 


40.7 


46.5 


52.1 


h(x!,y) 


o-+ 


11.8 


17.8 


23.5 


34.6 


34.6 


40.1 



Table 1: Mass spectrum in units of T for ++ and + glueballs corresponding to nor- 
malizable modes of fluctuations f(xi,y), g(xi,y), and h(xi,y). 

we get the relations 

9ttX ft 1 / N \ 

h(x 1: y). 



2 

9ym 



N ' 16tt 2 4 V2ttA, 
The supergravity fields f(xi,y) and h(xi,y) couple to the operators 



N 



2ir\J 5 



: TrF 2 , 



1 / N 



4 V2ttA 



TtFF. 



(2.21) 



(2.22) 



3 Mass Spectrum 



The glueball masses in QCD can be obtained by computing correlation functions of gauge 
invariant local operators or the Wilson loops, and looking for particle poles. Following 
[8, 9] correlation functions of local operators are related at large N and large gy M iV to 
tree level amplitudes of supergravity. We will consider the two point functions of the 
operators (2.22). The correlation function of O takes the form 



(O(x)O(y)) 



i 



-Mi\x-y\ 



(3.1) 



where O refers to (2.22), and Mi are the corresponding glueball masses. Therefore the 
spectrum of scalar glueballs 2 ++ and h can be obtained by finding the normalizable 
solutions to (2.11) for f(xi,y) and h(xi,y) respectively. The results are summarized 
in table 1, where the masses are given in units of T. We also included the spectrum 
associated with the field g(x 1: y). 

We can compare the ratio of masses of the lowest glueball states 0~ + and ++ with 
the lattice results [10] 

/Mq-+\ 

V Mq++ J supergravity 



= 1.20 



1.36 ±0.32, 



lattice 



(3.2) 



2 In the following we will use the notation J PC for the glueballs, where J is the glueball spin, and P, 
C refer to the parity and charge conjugation quantum numbers respectively. 
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and the results are in good agreement. 

The spectrum of masses related to g(xi,y) is higher than the spectrum of masses 
related to f(xi,y) since it corresponds to a higher dimensional operator. Unlike QCD 
where different operators with the same quantum numbers give the same mass spectrum, 
in the supergravity limit the corresponding mass spectrum is different. Since a TrF 4 
operator couples on the worldvolume theory to both / and g modes we expect that in the 
QCD limit the spectrum of masses associated with g will decouple in order for TrF 2 and 
TrF 4 to give the same spectrum. 

Let us make some comments on the results. In comparison to [4] the masses of 0~ + 
are lower. This is due to the fact that in [4] the computation of the spectrum of CT + was 
done using the RR 3-form which couples to an operator with 0~ + quantum numbers but 
of higher dimension than TrFF. An example of that is the ++ spectrum derived above 
using the / and g modes. Also, the masses of ++ are not in agreement with [4]. This 
is because the supergravity equation used in [4] to compute the ++ mass spectrum does 
not precisely correspond to TrF 2 which is the lowest dimension glueball operator with 
++ quantum numbers. 

Finally, in order to compute the coefficients q in (3.1) one has to evaluate the action for 
the supergravity fields requiring that they approach a boundary configuration at infinity. 
This is a procedure that can be done numerically, but will not be reported here. 



4 Topological Susceptibility 

The topological susceptibility Xt is defined by 

Xt = J d 4 x(TrFF(x)TrF F(0)). (4.1) 

The topological susceptibility measures the fluctuations of the topological charge of the 
vacuum. At large N the Witten-Veneziano formula [11, 12] relates the mass m v > to the 
topological susceptibility of Yang-Mills theory without quarks 

m 2 v , = —fxt- (4.2) 

J TV 

In this section we will compute the topological susceptibility in the supergravity de- 
scription. The discussion is closely related to the study of the 9 dependence of the vacuum 
energy in [13] . Working in Type IIA supergravity, the action of RR 1-form is given by 

I = ^T [ d w x^g-(d^C v - d u C,)(d,,C v , - d^C^g^'g™'. (4.3) 
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For the fluctuation C t = h(x\,y), the action simplifies to 

1 = A- / d'^^Ud^yW^K^v))^' 9 tt - (4.4) 

10 

The solution to the corresponding equation of motion (2.11) which has dependence only 
in the y coordinate is 

h(y) = h°°(l-%. (4.5) 

Plugging this into the action gives 

Using (2.6), (2.3), dropping the volume factor and differentiating twice with respect to 
h°° we get 

/ d^(d 4 (x)O 4 (0)) = U8X ^f T4 . (4.7) 
Recalling (2.22) and (4.1) we get the topological susceptibility 



5 Gluon condensate 



The gluon condensate (-j-z — F 2 ) plays an important role in the applications of Wilson 
OPE to the study of strong coupling phenomena in QCD via the condensate expansion. 
In this section we will compute the gluon condensate in the supergravity description. 

According to the supergravity/gauge theory correspondence (2.12), the one point func- 
tion of an operator corresponds to the first variation of the supergravity action. This 
quantity usually vanishes identically by the equation of motion. In asymptotically anti de 
Sitter spaces, however, there is a subtlety. The first variation is only required to vanish 
up to a total derivative term. Since asymptotically anti de Sitter space has a time-like 
boundary at infinity, one must carefully examine the possible boundary contribution. In- 
deed, the one point function of F 2 operator, corresponding to first variation of the action 
with respect to the field f(x±,y), will turn out to give a non- vanishing contribution. 

This result can be anticipated without any detailed computation. The non-extremal 
M5-branes has non-vanishing action density / (which is related to free energy F by I = 

64 

I = 2 6 3~ 7 n 3 N 3 V 5 T 5 = -n 3 N 3 V 5 T 5 . (5.1) 
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When the M5-brane is wrapped on a circle of period (A/A^)T~ 1 , we expect the corre- 
sponding free energy to be given simply by making the substitution V 5 — > V^T^X/N 
which leads to an expression for the free energy of non-extremal D4-branes: 

/ = 2 6 3-VAiVV 4 T 4 = ^7r 3 AiVV 4 T 4 . (5.2) 

The gluon condensate operator can be gotten by acting on the Born-Infeld action with 
A^. This suggests that one can act on / with the same operator to infer the value of 
(y4 — F 2 ) in the dual formulation: 

' F 2 } = = ^-n 3 XN 2 V 4 T 4 . (5.3) 



"Ag 2 M 1 d\ 2187 

6 Comparison 



We summarize the results obtained using the supergravity description: 

M ++ = 9.86T 



.3 



M -+ = 11. 8T 

167T 2 

QCD string tension = ~^~^ 2 
A 729 

1 647T 3 

TrF 2 (0)> = — — -iV 2 AT 4 . (6.1) 



^9ym K ~" 2187 
These results are the leading in the strong coupling expansion of the gauge theory. 

The lattice QCD results are [10] 4 



M 0++ ~ 1610 MeV 

M -+ ~ 2190 MeV 

QCD string tension ~ (440 MeV) 2 

Xt ^ (180 MeV) 4 

1 N 2 

(^^TrF 2 (0)> ~ ^-(400 MeV) 4 . (6.2) 

^9ym a 



3 For completeness we quote also the result for the QCD tension [15]. This result is corrected by a 
factor of 7r 3 compared to a previous version of this paper due to a numerical error in a previous version 
of [15]. 

4 The lattice results have systematic and statistical errors that we do not quote. 
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The entry for the gluon condensate is derived by scaling the N = 3 result of [16] by a 
factor of N/3. 

| .lGeV 4 = (^G 2 ) = (^F 2 ). (6.3) 

Note that plugging the phenomenological values Nf = 3, N = 3, ~ 1 GeV, f n ~ 
0.1 GeV in (4.2) leads to a prediction Xt ~ (180 MeV) 4 , in agreement with the lattice 
simulation. 

The supergravity results depend on two parameters, A and T. This should be thought 
of as the leading 1/A asymptotic behavior of the full string theory expression (F(X)T) d 
where d is the dimension of the observable. Ideally, the comparison with the lattice QCD 
result should be made in the limit where A is taken to be small and T is taken to be large. 
In this range of parameters, one can trust perturbative QCD, which tells us that F(\) 
goes as 

F(\) = Ce-V&. (6.4) 

where C is some numerical coefficient. The genuine QCD data is encoded in the ratio of 
this numerical coefficient C between different observables. Unfortunately, this is precisely 
the information that is not accessible via the supergravity approximation. In the absence 
of fully string theoretic formulation of this system that allows one to probe the small 
A region, one must contend oneself with qualitative comparisons. To facilitate such a 
comparison at a qualitative level, it is useful to go to the so-called "correspondence point" 
in the A-T parameter space. This is the point where dual descriptions of the theory 
crosses over from one to the other. As such, one can not trust neither description at a 
quantitative level, but one expects the results to match up to numerical factors of order 
one. One can estimate the cross-over point based on natural scale of our problem: A ~ ^ 
and T m Aqcd ~ 200 MeV. With this choice of A and T, the ratio of observables in units 
of energy comes out as follows: 



^STJGRA 

71 zf Lattice 
1VI 0++ 

M S _ U + GRA 

j^Lattice 



/ (QCD string tension) S UGRA \ 1/2 
\ (QCD string tension) 



Lattice 



/^UGRAXl/4 
I ,. Lattice I 

\Xt ) 
(Gluon condensate) SUGRA \ ^ 4 
, (Gluon condensate) Latticc ) 



= 1.2 

= 1.1 

= 2.0 

= 1.7 

= 0.9 
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It is encouraging to find that these ratio are indeed of order unity 5 . 

It is curious to note for comparison the known results in strong coupling lattice QCD 
[17]: 

M 0++ = -A\og(g 2 YM N)a- 1 
QCD string tension = \og(gy M N)a~ 2 . (6.5) 

They reflect the fact that the a' expansion is not the same as the lattice strong coupling 
expansion. In particular, while in the supergravity description the ratio of M ++ to the 
QCD string tension goes to zero, it diverges in the strong coupling lattice QCD. 

7 Summary and Discussion 

In this paper we computed various physical quantities which are of importance in QCD 
dynamics. The computations were performed in the framework of supergravity. Since 
the supergravity description is believed to provide the leading contribution in the strong 
coupling expansion of the gauge theory, our results should be viewed as a first step in 
understanding the QCD dynamics via the string/gauge theory correspondence. Clearly, 
a complete solution requires a detailed knowledge of the string dynamics with R-R back- 
ground. 

The dual supergravity description corresponds to taking the bare coupling to infinity. 
This is the opposite limit to continuum QCD where the bare coupling is taken to zero. 
Therefore it is not possible, for instance, to exhibit asymptotic freedom via the super- 
gravity description. In analogy with strong coupling lattice QCD, it is therefore far from 
clear what properties of continuum QCD should be captured by supergravity. The su- 
pergravity description includes also Kaluza-Klein modes that do not have corresponding 
degrees of freedom in QCD. They reflect the physics of the higher dimensions and should 
decouple in the small bare coupling limit. Since the background is singular in this limit, 
their decoupling is a subtle issue. 

We can repeat the analysis for QCD 3 . Here the procedure is to compactify on a circle 

the D3 branes of type IIB string theory and introduce anti periodic boundary conditions 

5 The numerical value for the ratio of the QCD string tensions is modified compared to a previous 
version of this paper in view of the change in the numerical value of the QCD string tension in (6.1). 
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for the fermions [2]. We find 6 



M ++ 



10.7T 



QCD string tension 



^9YM3 



1 



TrF 2 (0)> 




(7.1) 



The gluon condensate in QCD 3 has not been computed on the lattice yet. 

Acknowledgement 

We are grateful to O. Aharony, S. Carroll, C. Csaki, D. Gross, G. Horowitz, H. Ooguri, 
A. Pect, J. Polchinski, M. Srednicki, M. Teper, J. Terning, A. Tseytlin, and E. Witten 
for useful discussions. The work of AH was supported in part by the National Science 
Foundation under Grant No. PHY94-07194. 

References 

[1] J. Maldacena, "The Large N limit of superconformal field theories and 
supergravity," Adv. Theor. Math. Phys. 2 (1998) 231, hep-th/9711200. 

[2] E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in 
gauge theories," Adv. Theor. Math. Phys. 2 (1998) 505, hep-th/9803131. 

[3] D. J. Gross and H. Ooguri, "Aspects of large N gauge theory dynamics as seen by 
string theory," Phys. Rev. D58 (1998) 106002, hep-th/9805129. 

[4] C. Csaki, H. Ooguri, Y. Oz, and J. Terning, "Glueball mass spectrum from 
supergravity," JEEP 01 (1999) 017, hep-th/9806021. 

[5] M. J. Duff, H. Lu, and C. N. Pope, "The Black branes of M theory," Phys. Lett. 
B382 (1996) 73-80, hep-th/9604052. 

[6] S. S. Gubser, I. R. Klebanov, and A. A. Tseytlin, "Coupling constant dependence 

in the thermodynamics of N=4 supersymmetric Yang-Mills theory," Nucl. Phys. 

B534 (1998) 202, hep-th/9805156. 

6 For completeness we quote also the results for the ++ mass [4] and the QCD tension [15]. The 
spectrum in [4] was quoted in units of b, which in AdS„+i Schwarzschild space, is related to T by 
b = AttT/ti. 



12 



[7] I. Chepelev and A. A. Tseytlin, "Long distance interactions of D-brane bound 
states and longitudinal five-brane in M(atrix) theory," Phys. Rev. D56 (1997) 
3672-3685, hep-th/9704127. 

[8] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, "Gauge theory correlators from 
noncritical string theory," Phys. Lett. B428 (1998) 105, hep-th/9802109. 

[9] E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2 
(1998) 253, hep-th/9802150. 

[10] M. J. Teper, "Physics from the lattice: Glueballs in QCD: Topology: SU(N) for all 
N," hep-lat/9711011. 

[11] E. Witten, "Current Algebra Theorems for the U(l) Goldstone Boson," Nucl. Phys. 
B156 (1979) 269. 

[12] G. Veneziano, "17(1) Without Instantons," Nucl. Phys. B159 (1979) 213-224. 

[13] E. Witten, "Theta dependence in the large N limit of four-dimensional gauge 
theories," Phys. Rev. Lett. 81 (1998) 2862, hep-th/9807109. 

[14] S. W. Hawking and D. N. Page, "Thermodynamics of black holes in anti-de sitter 
space," Commun. Math. Phys. 87 (1983) 577. 

[15] A. Brandhuber, N. Itzhaki, J. Sonnenschein, and S. Yankielowicz, "Wilson loops, 
confinement, and phase transitions in large N gauge theories from super gravity," 
JEEP 06 (1998) 001, hep-th/9803263. 

[16] M. Campostrini, A. D. Giacomo, and Y. Gunduc, "Gluon Condensation in £77(3) 
Lattice Gauge Theory," Phys. Lett. B225 (1989) 393. 

[17] M. Creutz, Quarks, Gluons and Lattices. Cambridge University Press, Cambridge, 
UK, 1983. 



13 



